The next-to-leading-order QCD correction to inclusive J/tp{T) production in Z° decay 
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In this paper, we study the J/'i/'(T) production in Z boson decay in color-singlet model(CSM). 
We calculate the next-to-leading-order (NLO) QCD correction to Z — >■ Quarkonium + QQ, the 
dominant contribution in the CSM, with the vector and axial-vector parts in ZQQ vertex being 
treated separately. The results show that the vector and axial-vector parts have the same K factor 
(the ratio of NLO result to leading-order result) 1.13 with the renormalization scale fi=2mc and 
rric = 1.5GeV, and the K factor falls to 0.918 when applying the Brodsky, Lepage, and Macken- 
zie(BLM) renormalization scale scheme with obtained Hblm ~ 2.28GeV and rUc = l.SGeV. By 
including the contributions from the next-dominant ones, the photon and gluon fragmentation pro- 
cesses, the branching ratio for Z — >■ J/ipprompt -I- Jf is (7.3 ^ 10.0) x 10~^ with the uncertainty 
consideration for the renormalization scale and Charm quark mass. The results are about half of 
the central value of the experimental measurement 2.1xl0~*. Furthermore, the J/ip energy dis- 
tribution in our calculation is not well consistent with the experimental data. Therefore, even at 
QCD NLO, the contribution to Z — >■ J /ij^prompt + X from the CSM can not fully account for the 
experimental measurement. And there should be contributions from other mechanisms, such as the 
color-octet (COM) contributions. We define Rcs = '"rfglfj/'^x? ^^'^ obtain Rcc = 0.84 for only CSM 
contribution and Rcc ~ 0.49 for COM and CSM contributions together. Then Rcc measurement 
could be used to clarify the COM contributions. 

PACS numbers: 12.38.Bx, 13.38.Dg, 14.40. Pq, 12.39.Jh 



I. INTRODUCTION 

Heavy Quarkonium is an ideal system being used 
to study the perturbative and non-perturbative aspects 
of QCD. Firstly, the heavy quark mass sets a large 
scale for perturbative calculation. Secondly, the dilep- 
tonic decay of heavy quarkonium makes the identification 
and measurement efhcient. In 1995, the non-relativistic 
QCD(NRQCD), a rigorous effective theory in describing 
the production and decay of heavy quarkonium, was pro- 
posed [l|, and it makes the color-singlet model(CSM) [2| 
be its leading-order approximation in v (the velocity be- 
tween heavy quark and anti-quark in the meson rest 
frame) . More details on NRQCD and heavy quarkonium 
physics can be found in reference (3|. 

In recent years, there are many works on the next-to- 
leading-order(NLO) QCD correction for heavy quarko- 
nium productions. To explain the experimental measure- 
ment 4, 5] of JM production at the B factories, a series 
of calculations [a, |7| have been performed and revealed 
that the NLO QCD corrections can change the leading- 
order (LO) theoretical predictions considerably and the 
NLO results in CSM give the main contribution to the 
related processes. Together with the relativistic correc- 
tion 8^, it seems that all the experimental data for J/^p 
production at the B factories could be understood. For 
J/^ production in the hadron colliders, there are obvi- 
ously progress in the theoretical calculation. The NLO 
QCD correction to the CSM processes [9|, llO[ greatly en- 
hanced the pt (transverse momentum of J/ip) distribu- 
tion of J/iJj production at large pt region, and the pt dis- 



tribution of J/ip polarization is drastically changed from 
mostly transverse polarization at LO into mostly longi- 
tudinal polarization at NLO [10|. It is found that the 
NLO QCD correction to J/ip production for color-octet 
(COM) parts is quite small, about 10 percent 11|. Even 
including all these progresses, we still can not obtain an 
satisfactory explanation on both the pt distribution of 
the production and polarization for J/ip hadroproduc- 
tion. The partial next-to-next-to-leading-order(NNLO) 
calculations for T and J/ip hadroproduction show that 
the uncertainty form QCD higher order correction |12| is 
much bigger, but still can not cover the J/ip or T po- 
larization measurement. Recent studies reveal that the 
NLO QCD correction also plays an important role on 
J/ip pro duction at RHIC [l3J and the hadroproduction 
of Xc UM ■ The J/ip photoproduction once was considered 
as an positive example with the pt and z distribution well 
described by the NLO calculations in CSM ^J. But ei- 
ther the Pt distribution of the production or polarization 
for J/ip can not be well described by the recent NLO 
calculations in CSM [16[. It seems that the complete 
calculation at NLO in COM [17[ can account for the 
experimental measurements on the pt distribution. But 
the complete calculation on J/ip polarization at NLO in 
COM is a real challenge. 

With both the successful and unsuccessful aspects for 
theoretical progress in heavy quarkonium production, it 
is worthwhile to study more cases in detail. Such as J/ip 
production associated with photon [18|, QED contribu- 
tions in J/ip hadroproduction [19| , inclusive J/ip produc- 
tion from T decay J20| and J/ip production from Z de- 



cay. Heavy quarkonium production in Z decay has been 



widely studied in the CSM and COM at LO [2l|-|24|. and 
the measurement at the LEP by L3 Collaboration gives 
the branching ratio as [25| 

Br{Z -> J/lkprompt + X) 

= (2.1 ± 0.6(stai.) ± QA{sys.)tli{theo.)) x 10-^ (1) 



Br{Z -^ T(IS') + X) < AA x 10" 



(2) 



Theoretical investigation on this process indicates that 
even the dominant channel Z — >■ J /ij) + c + c in all the 
CSM ones at LO only gives the 1/3 prediction to the to- 
tal branching ratio of the experimental measurement. In- 
cluding the contribution of gluon fragmentation process 
in the COM can enhance the theoretical results about 
3 times [22|. This once is an evidence for the effect of 
the COM. It also have been studied in color-evaporation 
model(CEM) in reference [26| and obtained consistent re- 
sults with the experimental data. But the CEM always 
gives unpolarized J /tp in conflict with experimental mea- 
surements. By resuming the large logarithm from the 
large difference of J/ij] and Z mass, the COM prediction 
on J/i/; energy distribution((ir/rfz with z=2Ej/^/Mz) is 
roughly consist with the data [27| . Considering the larger 
impact of the NLO QCD corrections to the production of 
heavy quarkonium, it is necessary to investigate the NLO 
QCD correction to Z — > J/t/i ^ cc + X. In this paper, we 
calculate the NLO QCD correction to Z — >■ J/ip + cc + X, 
and also include the contributions from the gluon and 
photon fragmentation processes as well with only the 
CSM in consideration. The study could provide more in- 
sight to the effect of color-octet mechanism and put more 
constrains on the value of color-octet matrix elements. 

This paper is organized as follows. In Sec. II, we study 
the NLO QCD correction to the heavy quark association 
process with different schemes on the choice of renormal- 
ization scale. In Sec. Ill, we investigated dominant frag- 
mentation processes and give the total results on prompt 
J/ip production in Z decay. In Sec. IV, the summary 
and conclusion are presented. 




FIG. 1: Typical Feynman diagrams for J/^ + cc + X produc- 
tion in Z decay. B is the Born diagram. CV1~CV5 represent 
the counterterm diagrams and corresponding loop diagrams. 
V1~V7 represent the box and the anomalous triangle dia- 
grams. Rl and R2 are the real parts. 




FIG. 2: The typical amplitude squared diagrams for J/tjj + 
CC + X production in Z decay. 



II. THE HEAVY QUARK ASSOCIATION 
PROCESS 



For the calculation on Z — >■ J/tjj + cc + X at NLO, 
there are virtual and real correction parts as 



Z -^ J/ip + C + C 
Z ^ J/i} + c + c + g. 



(3) 

(4) 



There are vector and axial- vector parts in the coupling of 
Z boson to fermions, but the interference between them 
does not contribute in our calculation. Therefore, we 
study them separately. There are 4 Feynman diagrams 



for both parts at LO, 80 for the vector part and 76 for the 
axial- vector part at NLO. The typical diagrams are pre- 
sented in Fig. [TJ The dimensional regularization is used 
to regulate the the ultraviolet (UV) and infrared (IR) 
divergence, and the Coulomb singularity is regulated by 
introducing a small relative velocity between quark pair 
in the quarkonium and absorbed into the wave function 
of quarkonium. In calculating the axial-vector part, we 
have to face the 75 problem. The structure of all the 
amplitude squared diagrams could be classified into four 
cases shown in Fig. [5J 

Case 1. There are only one fermion-loop and two 75 
matrices appear in it. Then 75S can be moved together 



and give an identity matrix by 7I = 1. 

Case 2. There are two fermion-loops and the two 75 
matrices appear in one of them. It is the same as case 1. 

Case 3. There are two fermion-loops. Each of them 
has a 75. Because there are no UV and IR divergences 
in the loops, the dimension can be set as 4 safely. 

Case 4. The only special case are the two triangle 
anomalous diagrams. In this case we use the scheme de- 
scribed in reference [28| to handle it, which fixes the start- 
ing point to write down all the amplitude and abandon 
the cyclicity in calculating the trace of the fermion-loop 
with odd number of 75. This two triangle anomalous 
diagrams will not contribute at all according to Yang's 
theorem [29j when the two gluon lines are on mass-shell, 
but will contribute in our case since the two connected 
gluons are off mass-shell. 

The on-mass-shell (OS) scheme is used to define the 
renormalization constants Z^, -^2 and .Z3, which corre- 
spond to charm quark mass rric, charm field f/'c, and gluon 
field A^ while Zg for the QCD gauge coupling as is de- 
fined in the modified-minimal-subtraction(MS) scheme: 



wave functions at origin Rs , which is 
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(5) 



where /i is the renormalization scale, "fE is Euler's con- 
stant, /3o — ^Ca — ^TpUf is the one-loop coefficient 
of the QCD beta function and Uf is the number of ac- 
tive quark flavors. There are three massless light quarks 
u, d, s, and heavy quark c, so nf=A. In S'C/(3)c, color 
factors are given by Tp = ■^,Cf — ^,Ca = 3. And 
/3^ = /3o + iA/3)TF = (11/3)Ca - (4/3)rFn,/ where 
nif = Tiy — 1 = 3 is the number of light quarks flavors. 
Actually in the NLO total amplitude level, the terms 
proportion to SZ^ cancel each other, thus the result is 
independent of renormalization scheme of the gluon field. 
The above renormalization scheme and constant arc sim- 
ilar to those in reference [30[. The bottom quark should 
be considered for the calculation of T production. 

We use the Feynman Diagram Calculation(FDC) pack- 
age [3l| to generate Feynman diagram and amplitude, 
to do the tensor reduction and scalar integration, and 
to give the FORTRAN code for numerical calculation fi- 
nally. Because there are some large numbers generated in 
the program, the quadruple precision FORTRAN source 
is used. 

The leptonic width of J/ip{T) is used to extract their 



Tee = (1 - 



16q, 4a^e^ 



e{b) 



|i?f/^( 



T)|2 



J/^(T) 



J/^ 



Here the values of the parameters are chosen as Fee = 
5.55keV, FJ^ = 1.34keV [32[, a=l/137 and as = 
q,2;oop^277jq). The one- loop and two- loop running pro- 
gram of CTEQ6 are used to fix the LO and NLO values 
of as- The LO wave functions of heavy quarkonium are 
used to obtain the LO results in Fig. |31 IH E] and HI In 
the following calculation, a = 1/128 is used, and the cen- 
tral value of heavy quark mass is chosen as 77ic=1.5GeV 
(TO6=4.75GeV). We also use mc=1.4,1.6 GeV (TOh=4.65, 
4.85 GeV) for uncertainty estimate. The default choice 
of renormalization scale is 2me(2mi,) for J/ip (T). 

The LO and NLO partial decay width of Z -> J/ip + 
CC + X are presented in Table HI The difference between 
our LO results and the other LO results in literature is 
mainly due to the different choice of the wave functions 
at origin. The QCD correction enhance the partial de- 
cay width about 13% for both the vector part and the 
axial- vector part when the same wave function at ori- 
gin is used. This may provide a hint that the picture 
of heavy quark fragmentation into quarkonium works at 
these energy scale at NLO. It can also be seen that the K 
factors are insensitive to the variance of the quark mass. 
For the T production the similar results are presented in 
Table [III And it is easy to find that there is very small 
difference in K factors between the vector part and the 
axial- vector part. It could be thought as that the large 
bottom quark mass makes the fragmentation picture less 
effective than that in the J/ip production process. 

The renormalization scale dependence of the partial 
decay widths for J/ip and T are shown in Fig. |3] and 
Fig. m The QCD correction improve the scale depen- 
dence in small fi region and there are similar behavior 
for LO and NLO results in other region. In Fig. [5] 
and [6l the energy distribution of Jf-tjj and T are shown 
with z defined as 2i?j/^(x)A^z- The NLO QCD correc- 
tion shifts the maximum point of J/ip energy distribution 
from the large z region to the middle z region. But for 
T, the shifts is not manifest. 

To study the uncertainty from different choices of the 
renormalization scale, in addition to our default choice 
2me{mb) for J /'ip{T) in the calculation, we use other two 
schemes to fix the renormalization scale. At first, the 
decay width could be expressed as 



^NLO 



(A^) = r^^(A^)[i 



^iM(A + /3> ^ 



2mQ 



Bnf)].{6) 



Here the LO results depend on the renormalization scale 
through the running of the coupling constant. A and 
B are independent of the scale and /3q = 11 — 2nif/3. 
We extract the parameters in Eq. IH] and present them in 
Table inil 



m,(GeV) 


asifi) 


rl."'(keV) 


rW(keV) 


rW /r{0) 

^ V / ^ V 


T%(keY) 


rW(keV) 


r(i) /rC) 

'- AV/ ^ AV 


^ tot/ ^ tot 


1.4 


0.266 


19.6 


22.2 


1.13 


120 


136 


1.13 


1.13 


1.5 


0.259 


16.9 


19.1 


1.13 


103 


117 


1.13 


1.13 


1.6 


0.252 


14.8 


16.7 


1.13 


90.0 


102 


1.13 


1.13 



TABLE I: The partial decay width for J/tp with the renomialization scale /i = 2mc and different charm quark mass rric 



m6(GeV) 


as{^l) 


r<^'(keV) 


r«(keV) 


r«/rf 


rW(keV) 


r«(keV) 


p(i) /r(") 

^ AV/^ AV 


r(i)/r(0) 

-■- tot/ ^ tot 


4.65 


0.184 


5.50 


6.88 


1.24 


8.95 


11.1 


1.25 


1.24 


4.75 


0.183 


5.33 


6.68 


1.24 


8.61 


10.7 


1.25 


1.25 


4.85 


0.182 


5.17 


6.49 


1.24 


8.29 


10.3 


1.26 


1.25 



TABLE II: The partial decay width for T with the renormahzation scale fi = 2mb and different bottom quark mass mt. 




L<J / /-^Q 
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FIG. 3: The /i dependence of the partial decay width for 
Z — >■ J I'^cc + X with ntc = 1.5GeV and /io = 2mc. Here 
the LO results are calculated with the wave function at origin 
at LO and the a.a are fixed by one-loop running, and these 
choices are also applied for the LO plots in Fig.[4l [5]and[6l In 
all the figures, V presents the vector part result, A-V presents 
the ax;ial- vector part and total presents the sum of these two 
parts. 



m<:(GeV) 


r^/°^(keV) 


A 


B 


1.40 


176 


2.08 


-0.178 


1.50 


151 


2.12 


-0.182 


1.60 


131 


2.16 


-0.186 


mi,(GeV) 


r^°(keV) 






4.65 


17.8 


4.97 


-0.273 


4.75 


17.2 


5.05 


-0.275 


4.85 


16.6 


5.12 


-0.278 



TABLE III: The extracted parameters for Eq. [6] 



FIG. 4: The /i dependence of the partial decay width for 
Z -^ Tbb + X with mt = 4.75GeV and ^o = 2Mi,. 



Scheme I: From Fig. [3] and Fig. 21 it can be seen that 
there are the ^ points where the partial decay widths 
reach their maximum values. By using the Eq. [6l the 
values of ^ and partial decay widths can be obtained and 
presented in Table ITVl 

Scheme II: In Brodsky, Lepage, and Mackenzie (BLM) 
scheme [33|, the n// (light quark flavor) dependence of 
the QCD correction is absorbed into the running of as 
by shifting the renormahzation scale. An improved result 
on process e~^e~ -^ J/ipcc has been obtained in reference 
[7| . So we also try this scheme in our calculation and the 
results are presented in Table |V] and I VII It can be seen 
that the convergences of the perturbative expansions are 
all improved and the K factor is even lower than 1 for 
the J/ip production. 

The above two schemes give almost the same results 
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FIG. 5: The J/ij) energy distribution in Z — >■ J/^})cc with 
rric = 1.5GeV and ^ = 2Mc. 
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FIG. 6: The T energy distribution in Z — )■ Tbb with mt 
4.75GeV and ^ = 2Mt. 



m,(GeV) 


M(GeV) 


ryA°(keV) 


1.40 


2.26 


162 


1.50 


2.42 


139 


1.60 


2.58 


120 


m6(GeV) 




ry^°(keV) 


4.65 


6.48 


18.4 


4.75 


6.57 


17.8 


4.85 


6.66 


17.2 



TABLE IV: The maximum partial decay width for Z 
J/tp{T) + cc{bb) + X in the scheme I. 



m,(GeV) 


M*(GeV) 


Mf^*) 


r(o)(keV) 


r(i)(keV) 


p(i) /p(()) 


1.4 


2.14 


0.298 


176 


162 


0.919 


1.5 


2.28 


0.290 


151 


139 


0.918 


1.6 


2.42 


0.282 


131 


120 


0.918 



TABLE V: The partial decay widtli with different charm 
quark mass rric and renormalization scale /i = ^* in BLM 
scheme. 



m6(GeV) 


M*(GeV) 


CMsifJ.*) 


r(°'(keV) 


r(')(kev) 


p{i) /p(0) 


4.65 


6.18 


0.204 


17.8 


18.3 


1.03 


4.75 


6.29 


0.203 


17.2 


17.7 


1.03 


4.85 


6.39 


0.202 


16.6 


17.1 


1.03 



TABLE VI: The partial decay width with different bottom 
quark mass mt and renormahzation scale fi = /i* in BLM 
scheme. 



for both J/ip and T process. In the following discussion 
we will adopt the results from the BLM scheme. 



III. PHOTON AND GLUON FRAGMENTATION 
PROCESSES AND THE TOTAL RESULTS 

There are some QED processes which can give contri- 
butions comparable to that of the QCD ones in heavy 
quarkonium production 3J]. The contribution from the 
photon fragmentation processes was investigated in ref- 
erence [2J] and it gives non-ignorable contribution to the 
inclusive J/ip production in Z boson decay. Therefore, 
we further investigate the QCD correction to this photon 
fragmentation processes. At leading order, the following 
processes must be included. 



z -^ j/i' + r + r 

Z -^ J/ijj + q + q. 



(7) 
(8) 



Here l{q) is the lepton(quark) and the final results must 
be summed over e, fi and r(u, d, c, s, 6). We only pick out 
the photon fragmentation diagrams to calculate. These 
diagrams form a gauge invariance subgroup. All the typ- 
ical Feynman diagrams at LO and NLO are shown in 

Fig.m 

There are also the gluon fragmentation processes in 

CSM, 



Z -^ J/ip + q + q + g + i 



(9) 



Here the qq in the final states will be summed over 
u,d,c,s,b. Although they are at order aa^, the con- 
tribution of them is not too small [22, l23| . The typical 
Feynman diagrams are shown in Fig. [T] 

In evaluating these fragmentation processes, we set the 
renormalization scale as 2mc{2mi,). The NLO as and 




FIG. 7: The typical Feynman diagrams for the fragmentation 
processes on Z ^- J/ip + X. B is the Born diagram. CIl, 
Cql~Cq3 are the counterterm diagrams with it's correspond- 
ing loop diagrams for J/^ + l^l~ + X and J/ip + qq + X 
separately, V is the box diagram, R is the diagram for the 
real correction, and F is the gluon fragmentation process. 
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wave function for quarkonium are also used. Taking all 
the above processes in to account, we get the full results 
on the partial widths in Table I VIII and the energy dis- 
tribution in Fig. [5] 



m^(GeV) 


-pBLM 

J- J/i>+cc 


T^gluon 
^ QCD 


^ QED 


^ QED 


^ QED 


1.4 


162 


9.21 


10.5 


6.26 


4.36 


1.6 


120 


5.41 


8.12 


4.91 


3.43 



TABLE VII: The mass of charm quark is chosen as 1.4 GeV 
and 1.6 GeV, fi — ^blm for J /ijj -I- cc and /x = 2mc for other 
processes. Fq^*^ and Vqed present the contributions of the 
photon and gluon fragmentation processes respectively, (unit 
of decay widths: KeV) 



Combining all the above results together and timing a 
factor of 1.29 to include the contribution from the i/'' feed- 
down, we obtain the branching ratio of J /ijj production 
in Z decay as following: 



J/ip+cc+X 
NLO 


= (6.20 - 


^ 8.38) X 10~^ 


(10) 


^j.frag.pro. 


= (1.13- 


- 1.57) X 10"^ 


(11) 


Q^total 


= (7.33 r 


- 9.95) X 10-^ 


(12) 



Here we give the range of the branching ratio with the 
charm mass changing from 1.4 to 1.6 GeV. The total 
theoretical result is almost the half of the central value of 
the experimental measurement in Eq.(IT]). It is shown in 
Fig.[8]that the photon and gluon fragmentation processes 
contribute more in the lower energy region and the energy 
distribution can not fit the experimental data. 
Furthermore, we defined a ratio as 



Rcc — 



r{z ^j/ip + cc + x) 

T{Z ^ J/ij + X) ■ 



(13) 



Using the theoretical results obtained in the CSM, the 
ratio is about iJg."^*^ = 0.84 (0.85) with rUc = 1.4 (1.6) 



FIG. 8: The J/tp energy distribution in Z -;■ J/i/i -I- X with 
rric — 1.4 GeV, and jj. = hblm for J/ip + cc and n = 2mc for 
other processes. 



GeV. If we assume that the deriviation of the theoretical 
prediction from the central value of the experimental re- 
sults is from gluon fragmentation processes in the COM 
that was investigated in reference [2^ l23 1 , the ratio can 
be modified as 



p i 



TcsAiiZ ^ J/V^X)]},(14) 



where i?°g from gluon fragmentation processes in the 
COM is defined as 



no 



Y9- 



^r9^'Si(8)(^ ^ qq + j/^ + X)' 



(15) 



and the qq in the denominator are summed over u, d, 
c, s, b, and i?°g=0.17 is obtained from reference [23 1. 
Then we obtain i?;?/*^+^°*^ = 0.49 (0.41) for m^ = 
1.4 (1.6)GeV. The above analysis indicate that Rcc can 
be used to clarify the COM contribution. 



IV. SUMMARY AND CONCLUSION 

We have investigated all the processes that give main 
contributions to J/ip inclusive production in Z boson de- 
cay in the CSM. The results with NLO QCD correction 
are obtained. For the Z -^ J/^/j + cc process, the NLO re- 
sults only change the leading order results lightly, and the 
K factor is 1.13 with ii—2mc and insensitive to the charm 
quark mass. We also use two methods to estimate the de- 
pendence of the results on the choice of renormalization 
scale, and these two methods give almost the same partial 



decay width. The K factor even fall to 0.918 by using the 
BLM scheme. We also include the contributions of main 
fragmentation processes. The total branching ratio for 
Z ^ J/iP + X in CSM is (7.3 ~ 10.0) x IQ-s, about half 
of the central value of the experimental data 2.1x10^"*. 
We define Rcc — tiz*^j/1'x) ^^^ obtain R^c — 0.84 
for only CSM contribution and Rcc = 0.49 for COM 
and CSM contribution together. Then Rcc measurement 
could be used to clarify the COM contributions. In addi- 
tion the J /ip energy distribution is inconsistent with the 
experimental data too. But there are large uncertainties 
in the experiment results on the inclusive production of 
J /ip in Z decay, not only the total branching ratio but 
also the J /ip energy distribution. Further experimental 
measurement with more sample data is needed to clarify 
the situation. Maybe in the future Z factory these pro- 
cesses could obtain a detailed investigation. In the cal- 
culation, the K factor for vector and axial-vector parts 
of Z — > J /ip + CC + X are almost same. It may indicate 
that the mechanism of heavy quark fragmentation into 
quarkonium is dominant in this process even at NLO. 



ACKNOWLEDGMENTS 

We would like to thank Bin Gong and Hong-Fei 
Zhang for helpful discussion. This work is supported 
by the National Natural Science Foundation of China 
(No.10475083, 10979056 and 10935012) and by the Chi- 
nese Academy of Science under Project No. INFO-115- 
BOl, and by the China Postdoctoral Science foundation 
(20090460525). 



[1] G. T. Bodwin, E. Braaten and G. P. Lepage, Phys. Rev. 

D 51, 1125 (1995) [Erratum-ibid. D 55, 5853 (1997)] 

[arXiv:hep-ph/9407339 . 
[2] M. B. Einhorn and S. D. Ellis, Phys. Rev. D 12, 2007 

(1975); S. D. Ellis, M. B. Einhorn and C. Quigg, Phys. 

Rev. Lett. 36, 1263 (1976); C. H. Chang, Nucl. Phys. 

B 172, 425 (1980); E. L. Berger and D. L. Jones, Phys. 

Rev. D 23, 1521 (1981); R. Baier and R. Ruckl, Nucl. 

Phys. B 201, 1 (1982). 
[3] N. Brambilla et al. [Quarkonium Working Group], 

jarXiv:hep-ph/0412158; M. Kramer, Prog. Part. Nucl. 

Phys. 47, 141 (2001); J. P. Lansberg, Int. J. Mod. Phys. 

A 21, 3857 (2006). 
[4] K. Abe et al. [BELLE Collaboration], Phys. Rev. Lett. 

88, 052001 (2002); K. Abe et al. [Belle Collaboration], 

Phys. Rev. Lett. 89, 142001 (2002); K. Abe et al. 

[Belle Collaboration], Phys. Rev. D 70, 071102 (2004); 

P. Pakhlov et al. [Belle Collaboration], Phys. Rev. D 79, 

071101 (2009). 
[5] B. Aubert et al. [BABAR Collaboration], Phys. Rev. D 

72, 031101 (2005). 
[6] Y. J. Zhang, Y. j. Gao and K. T. Chao, Phys. Rev. Lett. 

96, 092001 (2006); Y. J. Zhang and K. T. Chao, Phys. 



Rev. Lett. 98, 092003 (2007); Y. J. Zhang, Y. Q. Ma and 
K. T. Chao, Phys. Rev. D 78, 054006 (2008); Y. Q. Ma, 
Y. J. Zhang and K. T. Chao, Phys. Rev. Lett. 102, 
162002 (2009); B. Gong and J. X. Wang, Phys. Rev. 
D 77, 054028 (2008); B. Gong and J. X. Wang, Phys. 
Rev. Lett. 100, 181803 (2008); B. Gong and J. X. Wang, 
Phys. Rev. Lett. 102, 162003 (2009); W. L. Sang and 
Y. Q. Chen, EXi^:0910.407ll [hep-ph]: D. Li, Z. G. He 
and K. T. Chao, Phys. Rev. D 80, 114014 (2009); 
Y. J. Zhang, Y. Q. Ma, K. Wang and K. T. Chao, Phys. 
Rev. D 81, 034015 (2010). 
[7] B. Gong and J. X. Wang, Phys. Rev. D 80, 054015 

(2009). 
[8] G. T. Bodwin, D. Kang, T. Kim, J. Lee and C. Yu, 
AIP Conf. Proc. 892, 315 (2007); Z. G. He, Y. Fan and 
K. T. Chao, Phys. Rev. D 75, 074011 (2007); G. T. Bod- 
win, J. Lee and C. Yu, Phys. Rev. D 77, 094018 (2008); 
Z. G. He, Y. Fan and K. T. Chao, Phys. Rev. D 81, 
054036 (2010); Y. Jia, arXiv:0912.5498 [hep-ph]. 
[9] J. M. Campbell, F. Maltoni and F. Tramontano, Phys. 
Rev. Lett. 98, 252002 (2007). 

[10] B. Gong and J. X. Wang, Phys. Rev. Lett. 100, 232001 
(2008); B. Gong and J. X. Wang, Phys. Rev. D 78, 
074011 (2008). 

[11] B. Gong, X. Q. Li and J. X. Wang, Phys. Lett. B 673, 
197 (2009). 

[12] P. Artoisenet, J. M. Campbell, J. P. Lansberg, F. Maltoni 
and F. Tramontano, Phys. Rev. Lett. 101, 152001 (2008). 

[13] S. J. Brodsky and J. P. Lansberg, Phys. Rev. D 81, 
051502 (2010); J. P. Lansberg, arXiv: 1003.4319 [hep-ph]. 

[14] Y. Q. Ma, K. Wang and K. T. Chao, arXiv: 1002.3987 
[hep-ph]. 

[15] M. 1. Kramer, Nucl. Phys. B 459, 3 (1996). 

[16] P. Artoisenet, J. M. Campbell, F. Maltoni and 
F. Tramontano, Phys. Rev. Lett. 102, 142001 (2009); 
C. H. Chang, R. Li and J. X. Wang, Phys. Rev. D 80, 
034020 (2009). 

[17] M. Butenschoen and B. A. Kniehl, Phys. Rev. Lett. 104, 
072001 (2010). 

[18] R. Li and J. X. Wang, Phys. Lett. B 672, 51 (2009); 
J. P. Lansberg, Phys. Lett. B 679, 340 (2009). 

[19] Z. G. He, R. Li and J. X. Wang, arXiv:0904.1477 [hep- 
ph]; Z. G. He, R. Li and J. X. Wang, Phys. Rev. D 79, 
094003 (2009). 

[20] Z. G. He and J. X. Wang, Phys. Rev. D 81, 054030 
(2010). 

[21] B. Guberina, J. H. Kuhn, R. D. Peccei and R. Ruckl, 
Nucl. Phys. B 174, 317 (1980); W. Y. Keung, Phys. 
Rev. D 23, 2072 (1981); K. J. Abraham, Z. Phys. C 44, 
467 (1989); V. D. Barger, K. m. Cheung and W. Y. Ke- 
ung, Phys. Rev. D 41, 1541 (1990); K. Hagiwara, 
A. D. Martin and W. J. Stirling, Phys. Lett. B 267, 527 
(1991) [Erratum-ibid. B 316, 631 (1993)]; E. Braaten, 
K. m. Cheung and T. C. Yuan, Phys. Rev. D 48, 
4230 (1993); J. Jalihan-Marian, arXiv:hep-ph/9401229 
P. Ernstrom, L. Lonnblad and M. Vanttinen, Z. Phys. C 
76, 515 (1997); G. A. Schuler, Int. J. Mod. Phys. A 12, 
3951 (1997). 

[22] K. m. Cheung, W. Y. Keung and T. C. Yuan, Phys. Rev. 
Lett. 76, 877 (1996); P. L. Cho, Phys. Lett. B 368, 171 
(1996); 

[23] S. Baek, P. Ko, J. Lee and H. S. Song, Phys. Lett. B 389, 
609 (1996); 

[24] S. Fleming, Phys. Rev. D 48, 1914 (1993). 



[25] M. Acciarri ei al. [L3 Collaboration], Phys. Lett. B 453, 

94 (1999). 
[26] E. M. Gregores, F. Halzen and O. J. P. Eboli, Phys. Lett. 

B 395, 113 (1997). 
[27] C. G. Boyd, A. K. Leibovich and I. Z. Rothstein, Phys. 

Rev. D 59, 054016 (1999). 
[28] J. G. Korner, D. Krcimer and K. Schilcher, Z. Phys. C 

54, 503 (1992). 
[29] C. N. Yang, Phys. Rev. 77, 242 (1950). 
[30] B. Gong and J. X. Wang, Phys. Rev. D 77, 054028 



(2008). 
[31] J. X. Wang, Nucl. Instrum. Meth. A 534, 241 (2004). 
[32] C. Amsler et al. [Particle Data Group], Phys. Lett. B 

667, 1 (2008). 
[33] S. J. Brodsky, G. P. Lepage and P. B. Mackenzie, Phys. 

Rev. D 28, 228 (1983). 
[34] K. Y. Liu, Z. G. He and K. T. Chao, Phys. Rev. D 68, 

031501 (2003). 



